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ABSTRACT 

We report the detection of spatially resolved C0(J=1— >-0) emission in the z~3.4 submillimeter 
galaxies (SMGs) SMM J09431+4700 and SMM J13120+4242, using the Expanded Very Large Array 
(EVLA). SMM J0943 1+4700 is resolved into the two previously reported millimeter sources H6 and 
H7, separated by ~30kpc in projection. We derive C0(J=1— )-0) line luminosities of ^^0(1-0) ~ 



(2.49 ±0.86) and (5.82 ± 1.22) x lO^o K kms"i pc^ for H6 and H7, and ico(i-o) = (23.4 ±4.1) x 

10i°K kms^ipc^ for SMM J13120±4242. These are '-l.5-4.5x higher than what is expected from 
simple excitation modeling of higher- J CO lines, suggesting the presence of copious amounts of low- 
excitation gas. This is supported by the finding that the CO(J=1-^0) hne in SMM J13120+4242, 
the system with lowest CO excitation, appears to have a broader profile and more extended spatial 
structure than seen in higher- J CO lines (which is less prominently seen in SMM J09431±4700). 
Based on ico(i-o) 



and excitation modeling, we find Mgas=2. 0-4.3 and 4.7-12.7 x 10^" for H6 and 



H7, and Mgas=18.7-69.4x 10 



10 



Mq for SMM J13120±4242. The observed CO(J=1-^0) properties are 
consistent with the picture that SMM J09431±4700 represents an early-stage, gas-rich major merger, 
and that SMM J13120+4242 represents such a system in an advanced stage. This study thus highlights 
the importance of spatially and dynamically resolved C0(J=1— ?>0) observations of SMGs to further 
understand the gas physics that drive star formation in these distant galaxies, which becomes possible 
only now that the EVLA rises to its full capabilities. 

Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

With star formation rates of typically >500Moyr~^, 
submillimeter galaxies (SMGs) are the most intense star- 
burst s know n at early cosmic times (see review by Blain 
et al.|2002)- A significant fraction of the stellar mass as- 
sembly in massive galaxies at high redshift is thought 
to occur at early epochs (z>2) on rapid timescales 
(<100Myr) throu gh a dusty, gas-rich SMC phase (e.g., 
Greve et al. |2005() . making SMGs the likely progenitors 
of massive galaxies in the present day universe. 

A particularly useful probe to better understand 
the physical conditions for star formation in these ex- 
treme starbursts are studies of the molecular interstellar 
medium that fuels the formation of stars. Dense molec- 
ular gas traced through CO was detected in 34 SMGs 
to date, typically revealing gas reservoirs of few 10^° Mq 
(see review by Solomon & Vanden BoutL2005). However, 
most of these detections are of CO rotational lines from 
J>3 transitions, which may bias these studies toward 
highly-excited gas, and do not necessarily trace the en- 
tire molecular gas reservoir as traced by C0(J=1— ?>0). 
Indeed, recent studies hint at the presence of substantial 
amounts of low-excitation gas i n som e SMGs (e.g., Hain- 
line et al. ,200& : Carilli et al. IMl Harris et al. I20T0I: 
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Ivison et al. I2010t Riechers et al. I2010a[) . Unfortunately, 
previous C0(J=1— >-0) studies in SMGs (which do not 
suffer from excitation bias) were either carried out with 
single-dish telesopes lacking spatial information, or with 
the old correlator of the Very Large Array (VLA) lacking 
detailed spectral information. 

To overcome the limitations of previous studies, we 
have initiated a systematic study of the CO(J=1^0) 
content of SMGs and other galaxy populations w ith th e 
Expanded Very Large Array (EVLA; Perley et al. 120111) . 
We here report early results on two z~3.4 SMGs that 
appear to trace different stages of gas-rich, major mergers 
only -1.9 Gyr after the Big Bang. SMM J09431+4700 
(z=3.346) hosts two SMGs detected in J>4 CO emission, 
separated by ~30kpc in projection, likely repres enting 
an early phase of a merger (e.g., Tacconi et al. 120061: 
Engel et al. lMol) . SMM J13120+4242 (z=3.408) shows a 
complex, disturbed structure in C0( J=6— ?>5), suggesting 
that it may represent an advanced-stage merger (Engel et 
al. [20101). We use a concordance, flat ACDM cosmology 
throughout, with gn= 71 knis~^ M pc^-^, r2M=0.27, and 
r2A=0.73 (Spergel et al. [200l [2001 . 

2. OBSERVATIONS 

We observed the CO(J=1^0) {vrcst = 115.2712 GHz) 
emission hne toward SMM J09431+4700 and 
SMMJ13120+4242, using the EVLA. At z=3.346 
and 3.408, this line is redshifted to 26.5235 and 
26.1505GHz ('-1.14cm; K band). Observations were 
carried out under good weather conditions in D array 
on 2010 May 19 and in C array between 2010 December 
16 and 2011 January 19, resulting in 4.0 hr (2.0 hr) 
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Fig. 1.— EVLA maps (left) and spectra (right) of CO(J=1^0) emission toward the z=3.346 SMG J09431+4700 H6 (top) and H7 
(bottom). The resolution of both maps is 0.91" xO. 81". The crosses indicate continuum positions reported by Tacconi et al. (2006). 
Contours are shown in steps of la=30 /ijybeam"^, starting at ±2(t. The spectra (histograms) are shown at 45 kms~^ (4 MHz) resolution. 
Zero velocity corresponds to z=3.346. The curves indicate Gaussian fits to the spectra. Top: Emission toward H6, averaged over 475 kms^^ 
(42 MHz; central velocity: 585 kms~^). Bottom: Emission toward H7, averaged over 543 kms~^ (48 MHz; central velocity: —123 kms^^). 



total (on-source) observing time for SMM J0943 1+4700 
(C array), and 7.5 hr (4.6 hr) for SMMJ13120+4242 
(C+D array). The nearby quasars J0920+4441 and 
J1312+4828 were observed every 4 to 7 minutes for 
pointing, secondary amplitude and phase calibration. 
For primary flux calibration, the standard calibrator 
3C286 was observed, leading to a calibration that 
is accurate within '^10%. Observations were set up 
using a total bandwidth of 252 MHz (dual polarization; 
corresponding to ^2900 kms^^ at 1.14cm) with the 
WIDAR correlator. 

For data reduction and analysis, the AIPS package was 
used. All data were mapped using 'natural' weighting. 
The data result in a final rms of 105 and 70 /iJy beam~^ 
per ~45 kms'^ (4 MHz) channel for SMM J09431+4700 
and SMM J13120+4242, respectively. Maps of the 
velocity-integrated CO ,7=1^-0 line emission yield 
synthesized clean beam sizes of 0.91" x 0.81" and 
1.01"x0.72" at rms noise levels of 30 and 12/.tJy beam~^ 
over 543 and 1513 kms^^ (48 and 132MHz) for 
SMM J09431+4700 and SMM J13120+4242. 



3. RESULTS 



3.1. 



i700 



SMM J 09431 - 

We detected spatially resolved CO(J=1^0) emission 
toward SMM J09431+4700. The emission is resolved into 
two components that are separated by ~30kpc in projec- 
tion, and are identified with the radio/millimeter contin- 
uum sources H6 and H7 (Fig.[T]). Component H7 appears 
marginally resolved at the '^6.5kpc resolution of our 
observations, suggesting an approximate source radius 
of -3±lkpc. This suggests that the CO(J=1->0) line 
emission may be somewhat more extended than what is 
observed in higher-J lines (1.85±0.60 and 1.13±0.23kpc 
in the CO J=4-S>3 and 6-)^5 lines; Engel et al. |201C1[) . 
Component H6 appears unresolved at the moderate 
achieved signal-to-noise ratio. There is marginal evi- 
dence for some low surface brightness emission in be- 
tween H6 and H7, but more sensitive observations are 
required to investigate this in more detail. We do not de- 
tect continuum emission toward either component down 
to a 3(7 limit of 44/iJybeam~^ at 1.1cm (rest- frame 
2.6 mm), consistent with the spectral energy distribution 
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Fig. 2.— EVLA maps (left) and spectra (right) of CO( J=l^-0) emission toward the 2=3.408 SMG J13120+4242. The r esolut ion of both 
maps is 1.01" xO. 72". The cross indicates the center of the CO(J=4— >3) emission detected at ^6" resolution (Greve et al. 120051 1. Contours 
are shown in steps of lcr=12 and 29 /xjy beam~^ (top/bottom), starting at ±2cr. The spectra (histograms) are shown at 46 kms~^ (4MHz) 
resolution, and are extracted over the emission regions delimited by the 2cr contours in the maps to the left. Zero velocity corresponds to 
2=3.408. The solid curves indicate Gaussian fits to the spectra. Top; Total emission over 1513 kms~i (132 MHz). The empty histo gram 
and dashed curve in the right panel indicate the single-dish measurement with the Green Bank Telescope at ~28" resolution (Hainline et 
al. 120061 ). Bottom: Emission over the central 344 kms ~^ (30 MHz). The empty histogram shows the CO(J=4— )-3) emission line, scaled 
down by a factor of 5 in flux density (Greve et al. 120051 ). 

of this source. 

From Gaussian fitting to the hne profiles of 
H6 and H7, we obtain line peak strengths of 
5'i.=170±58 and 232±47AiJy at line FWHMs of 
dii=339±136 and 580±143 kms^^, centered at red- 
shifts of z=3.3545±0.0009 and z=3.3442±0.0009, re- 
spectively. The line widths correspond to 1.9 x 
(H6) and 1.5x (H7) those measured in CO(J=6^5) 
(Engel et al. |2010( ). The line parameters for 
H6 and H7 correspond to velocity-integrated emis- 



sion line strengths of /( 
0.143±0.030Jy kms 
minosities of L'r 



CO(l-O)- 



=0.061±0.021 and 



suggesting CO(J= 
and 



4-^0) line lu- 



co-(2.49±0.86) and (5.82±1.22) x 
10^° (/iL/1.2)"^K kms"^pc^ (where /LtL=1.2 is the lens- 
ing magnification factor; this SM G is lo cated behind the 
galaxy cluster A851; Cowie et al. I2002f ). 

3.2. SMMJI312O+4242 

We detected spatially resolved C0(J=1— J'O) emission 
toward SMM J13120-K4242. The emission peaks on a 



(2.0"±0.6")x(0.9"±0.4") (15.0x6.7kpc2) size region 
(Fig. [2 bottom). The line profile averaged over this 
region is consistent with that seen in the C0(J=4— s>3) 
line (empty histogram in Fig. [31 bottom; Greve et al. 
120051 ). and its complex morphological structure is con- 
sistent with what is seen in C0(J=6— >5) line emission 
(Engel et al. 120101 ). This component, however, carries 
only <40% of the observed CO(J=1->0) line flux, and 
appears to be embedded in a kinematically a few times 
broader component that extends out to >2x larger 
scales (Fig. O top). The line profile over this whole 
extended structure appears to consist of multiple peaks. 
For consistency, we fit the line with a Gaussian profile, 
yielding S^=580±65 fiJy and d?;=1153±235 kms^^, 
and /co(i-o)=0-464±0.081Jy kms~^, centered at 
z=3.4078±0.0014. Within the relative uncertainties, 
the line profile and parameters are consistent with a 
single-dish measurement of this line, carried out with 
the Green Bank Telescope at 28" resolution (empty 
histogram and dashed curve in Fig. [2l top; Hainline et 
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Fig. 3.— CO excitation ladders (points) and LVG models (lines) for SMM J13120+4242 and SMM 09431+ 4700 H6 and H7. The 
CO(J=4^-3) and CO(J=6-*-5) data are adopted from Grcvc ct al. (2005), Tacconi et al. (2006) and Engel et al. (|20TOl ). The insets show 
zoomed-in versions of the model fits close to the CO(J=l— >0) line. The data for all sources are fit well by two-component LVG models, 
as indicated by the solid lines. The dashed lines indicate the two components of each model, represented by a low-excitation component 
with a kinetic temperature of Ti;in=25K and a gas density of Pgas = 10^ cm~'^, and a high excitation component with Ti5i,i=40K and a gas 
density of pgaa=10*'^ cm""^ (this component has a 1.26x lower pgaa for H7). The high-excitation components contribute 20%, 65%, and 
49% to the observed CO(J=l— >0) intensities. For comparison, the dotted line in the left panel shows a single-component fit to the J>1 
hues with Tkin=40K and pgas=10'*" cm^^. 

Pgas=10^'^ cm""^, and a more 'dense', high-excitation 
component with Tkin=40 K (comparable to th e dust tem- 
perature in SMM J13120-h4242; Kovacs et al. |2006[) and 



al. I2006D . This suggests that the large apparent spatial 
extent of the CO(J=1^>0) emission is real, and that 
virtually all line flux is recovered in the EVLA map. 
The line FWHM corresponds to 2.2 x and 1.3 x those 
measured in the C0(J=4— >-3) and C0(J=6— J^S) lines. 
We do not detect any rest-frame 2.6 mm continuum 
emission down to a Scr limit of 60 /iJy beam~^, consistent 
with the spectral energy distribution of this source. 
Averaged over the full region where C0(J=1— t'O) line 
emission is detected, this suggests a conservative 3ct 
limit of ~270 fiJy. From the line intensity, we derive 
L'pQ=(23.4±4.1) X 10i°K kms-ipc^. 

4. ANALYSIS 
4.1. Line Excitation Modeling 

Besides the EVLA CO(J=1^>0) observations pre- 
sented here, SMM J09431-H4700 and SMM J13120-f 4242 
were studied in C0( J =4— >-3) and C0( J=6^>5) line emis- 
sion (see Engel et al. 120101 for a summary of previous 
studies). Based on the observed CO excitation ladders, 
we can constrain the line radiative transfer through Large 
Velocity Gradient (LVG) models, treating the gas ki- 
netic temperature and density as free parameters. For 
all calculations, the H2 ortho-to-para ratio was fixed 
to 3:1, the cosmic microwave background temperature 
was fixed to 11.84 and 12.01 K (at z=3.346 and 3.408), 
and the Flower (2001) CO collision rates were used. 
We adopted a CO abundance per velocity gradient of 
[C01/ ( d^/dr ) = 1 X 10-5pc (kms-i)-i (e.g., Weifi etaL 
I2005al 120071: Riechers etaL I2006D . Components H6 and 
H7 of SMM J09431-h4700 are fitted individually. 

The sources are poorly fit by single-component mod- 
els, which underpredict the observed C0( J=l— ?>0) fluxes 
by factors of ~ 1.5-4. 5 when fltted to the higher- J 
linesO The data for all sources can be fltted reason- 
ably well with two gas components, which are repre- 
sented by a 'diffuse', low-excitation component with a 
kinetic temperature of rkin=25K and a gas density of 



Pgas=10 cm"-^ in our models (pgas=10 cm~'^ is used 
for H7, as it provides a slightly better flt). The two 
gas components differ in relative strength between the 
sources (Fig. For H6 and H7, we flnd that the 

dense gas components have surface flUing factors of 13% 
and 7% relative to the low-excitation component. For 
SMM J13120-h4242, we flnd that only -2% of the surface 
area associated with emission from the low-excitation 
component also shows emission from the dense gas com- 
ponent. This suggests that only a small fraction of the 
volume in these galaxies occupied by gas is associated 
with dense regions. Also, these model parameters sug- 
gest that the low-excitation gas components flll a domi- 
nant fraction of the area over which the sources are re- 
solved in these observations. 

The physical properties of the low-excitation gas com- 
ponents in these models are comparable to those of the 
gas in nearby spiral galaxies and 'normal' high-z star- 
forming galaxies (e.g., Dannerbauer etal. 2009) and the 
low-excitatio n com ponents found in other SMGs (e.g., 
Carilli et al. 120101 : Riechers et al. 2010a). The proper- 
ties of the high-excitation components are also compara- 
ble to what is found in other SMGs, ultra-luminous in- 
frared galaxy (ULIRG) nuclei, and hi gh-z FIR-luminous 
quasars (e.g., Riechers et al. 120061 I2009i Weifi et al. 
2005b). However, it appears that the relative fractions 
of low- and high-excitation gas differ between individual 
systems, perhaps due to different evolutionary stages. 

4.2. Gas Masses and Surface Densities, Dynamical 
Masses, Gas Mass Fractions 

Based on the observed C0(J=1— ?>0) line luminosities, 
we can derive the total masses of the gas reservoirs 
in our targets. For SMGs, a ULIRG-like conversion 
factor of aco=A^gas/ico=O-8^0 (K kms 



IS 



^ Given the limited constraints on H6, 
marginally fit with a single component. 



commonly adopted (e.g., Downes & Solomon 119981 : 
this source can be Tacconi et al. I2008D . This suggests Mgas=2.0 and 
4.7x1010 (/LtL/1.2)-iM0 for SMM J0943 1-^4700 H6 and 
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Fig. 4.— Overlays of the CO{ J=1^0) emission (contours, as in Figures[T]and[2ll toward J09431+4700 (left) and the center of J13120+4242 
(right) on Hubble Space Telescope Advanced Camera for Surveys (ACS) F 814W (McGrath et al. 2008) and F775W images (gray scale; 
rest-frame ~180nm; smoothed with a O.l" Gaussian kernel; Hainline et al. 12006 1. The gray and black contours in the left panel indicate 
the velocity ranges for H6 and H7 as in Fig. [T] 



H7, and Mgas=18.7x 10^" M© for SMM J13120+4242, 
respectively. However, a low, ULIRG-like conversion 
factor may not be appropriate for the low-excitation 
gas components suggested by our models. Assum- 
ing aco=3-5M0 (K kms~^ pc^)"-'^ for these compo- 
nents (as suggested for gas components with com- 
parable excitation fou nd in disk-like high-z galax- 
ies; e.g., Daddi et al. I2010D . we find Mgas=4.3 and 
12.7X IQi" (AiL/1-2)-^ Mq for SMM J09431+4700 H6 and 
H7, and Mgas=69.4x 10^° Mq for SMM J13120+4242, re- 
spectively. 

Assuming FWHM radii of 3 kpc for SMM J0943 1-^4700 
H6 and H7, the gas masses correspond to average gas sur- 
face densities of i;gas=700-1500 and 1600-4500 Mq pc^2_ 
The CO{J=l^Q) emission in SMM J13120-H4242 is dis- 
tributed over a projected area of ^200 kpc^, correspond- 
ing to Egas=96O~36OOM0pc-2. 

SMM J13120+4242 has a stel lar mass of 
13.5xl0i"Afo (Engel et al.|2010D. This yields a baryonic 
gas mass fraction of /^a"^=Mgas/(Mgas+Af*)=0. 58-0.87, 
where the spread of values represents the difference in 
conversion factors as outlined above. 

Usi ng the "isotropic vi rial es timator" (e.g., Tacconi 
et al. I200I Engel et al. [20Tol ). we can also estimate 
the dynamical masses of our targets. Assuming FWHM 
radii of 3 kpc for SMM J0943 1-^4700 H6 and H7, we find 
Mdy„=9.7 and 28.3xlOi°M0, respectively. This corre- 



sponds to gas mass fractions /gas=Afgas/Mdyn=0.21-0.45 
and 0.16-0.45 for H6 and H7, where the spread of val- 
ues represents the use of different conversion factors as 
outlined above. This lies within the typical ra nge of 
value s estimated for SMGs (e.g., Tacconi et al. I2006L 
[2008[) . The CO(J=1^0) emission in SMM J13120-F4242 
is spread over an area with a characteristic radius of 
7.9 kpc. Assuming that the multiple peaks can be de- 
scribed by a single potential that is represented correctly 
by the FWHM on the Gaussian fit used above, we find 
Mdyn=29OxlOi°M0 for SMM J13120-h4242. This cor- 
responds to /gas=0.06-0.24. Under the assumption that 
the system is gravitationally bound, and that Afgas+Af^, 
corresponds to the bulk of baryonic matter in this sys- 
tem (the contributions of dust and black hole mass are 
expected to be minor) , this would imply that this SMG is 



dominated by dark matter. However, given the complex 
velocity structure of this system, we consider it likely 
that A/dyn is overestimated by a factor of 2-3 with this 
simplified approach. Thus, we consider the gas fraction 
for SMM J13120-1-4242 a lower limit. 

4.3. Origin 0/ the Molecular Line Emission 

SMM J13120-t-4242 and both components of 
SMM J09431+4700 are associated with galaxies de- 
tected with the Hubble Space Telescope in rest-frame 
optical/ultraviolet light (Fig. |4]). H7, the more CO- 
luminous component of SMM J13120-h4242, is not 
significantly d etecte d in i?-band images of this source 
(Tacconi et al. l2Q06D . which may indicate that it is more 
obscured than H6. The CO emission centroids in H6 and 
H7 are offset by ■~1.5kpc (0.2") from the centers of the 
optical emission, but are clearly associated with fainter, 
potentially more obscured regions in the galaxies. The 
optical emission in SMM J13120-I-4242 appears to be 
associated with one of the CO peaks in the central 
region, and thus, potentially with one component of a 
major merger (or a region of low dust obscuration). 

5. DISCUSSION 

We have detected spatially resolved C0( J=l— s^O) emis- 
sion toward the z-3.4 SMGs SMM J09431-f-4700 and 
SMM J13120-I-4242. We resolve the gas reservoirs down 
to <5kpc scales (given the sizes of our synthesized 
beams). CO(J=1^0) emission in SMM J09431-h4700 
is detected toward both radio/millimeter continuum 
sources H6 and H7, which have been detected in higher- 
J CO lines in earlier work. There is marginal evidence 
for additional C0(J=1— )>0) emission in between the two 
sources, which are separated by ~30kpc in projection. 
This, and their proximity in redshift (dz<0.01) is con- 
sistent with the picture that this SMG is undergoing an 
ear ly sta ge massive, gas-rich merger (see also Engel et 
aL [20T0l ). The CO(J=1-^0) line in H7, the more CO- 
luminous component, appears to show a ~50% larger line 
width and size relative to what is measured in higher- J 
CO lines, but the values are marginally consistent within 
the relative uncertainties. 

The CO(J=1^0) emission in SMM J13120-1-4242 
shows a bright central peak that appears to be associated 
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with the emission that dominates the flux in higher- J 
CO hnes. It also shows >2x broader emission (as also 
seen in single- dish o bservations at low spatial resolution; 
Hainline et al. 120061 ) distributed over several peaks that 
extends out to >15kpc scales, but appears to be faint in 
higher- J CO hnes. Thus, the EVLA C0( J=l->0) obser- 
vations reveal a spatially resolved CO excitation gradient 
in this SMC This may suggest that this SMC is in an 
advanced stage of a merger, where the warm and dense, 
highly excited gas component predominantly arises from 
a starburst in the "overlap" region of two gas-rich galax- 
ies, and the broader, more extended, less highly excited 
gas is associated with tidal structure redistributed by 
the merger. More sensitive observations are required to 
study the dynamical structure of this complex system in 
more detail. 

The CO(J=1-J-0) properties of SMM J13120-f 4242 
contrast those of high-z quasars, which appear to be 
dominated by less extended (~5kpc), highly excited 
molecular gas reservoirs, with only little evide n ce for 
low-excitation gas (e.g^ Riechers et al. 120061 : 120091 : 
I2011al: Weifi et al. I2007D . They also contrast those of 
lensed z~3 Lyman-break galaxies, which typically have 
an order of magnitude less massive, com pact (^ 1-2 kpc) 
C0(J=1— >0) reservoirs (Riechers et al. I2010bl: see also 
highe r- J C O observations by Baker et al. 120041 : Coppin 
et al. I2007D . Massive, gas-rich star- forming galaxies can 
host gas reservoirs comparable in extent to those found 
in our target SMGs, but their spatial and velocity struc- 
ture typicall y app ears more symmetric and ordered (e.g., 
Daddi et al. I2O1O0 . Also, their overa ll CO excitation ap- 
pears to be lower (Dannerbauer et al. 120091 : Aravena et al. 
|201C1[) . Low CO excitation was also found in the z—lAA 
SMC HRIO, which was initially selected as an Extremely 
Red Obj ect (E RO: Andrea ni et al. [2000; Papadopoulos 
& Ivison[200| Greve et al. l200l . 

The new CO(J=1-^0) observations of the z-3.4 SMGs 
SMM J09431+4700 and SMM J13120+4242 provide sup- 
porting evidence that they represent early and advanced 
stage, gas-rich major mergers. Recent, lower resolution 



C0(J=1— J>0) studies of the ~20kpc separ ation merger 
SMM J123707+6214 (z=2.488; Ivison et aL [20Tll: Riech- 
ers et al. l2011bD may suggest an intermediate merger 
stage (with two separated components, but closer in pro- 
jected distance and redshift than in SMM J09431-I-4700), 
in between those observed in these two targets. These 
EVLA observations thus may be a first step toward es- 
tablishing a molecular gas-based "merger sequence" for 
gas-rich starburst galaxies at high redshift, which may 
crucially constrain formation models of SMGs (e.g., Dave 
et al. 2010; Hayward et al. 2011). 

If the C0(J=1— >0) properties of the 2;~3.4 systems 
studied here are representative of their parent popula- 
tion, SMGs in advanced merger stages exhibit substan- 
tially more extended, broader low- J CO line emission 
than in earlier merger stages. Observations of high-J 
CO line emission alone will underestimate the gas con- 
tent and extent of the tidal structure in advanced-stage 
mergers, yielding an incomplete picture of the processes 
that drive the gas dynamics in such systems. 

The observations presented here thus show the 
key importance of spatially and dynamically resolved 
CO(J=1^0) observations of SMGs to understand the 
gas physics that drive star formation in these luminous, 
massive gas-rich high redshift galaxies. In combination 
with higher- J CO line observations with shorter wave- 
length interferometers such as the upcoming Atacama 
Large (sub) Millimeter Array (ALMA), this makes the 
EVLA a uniquely powerful tool to distinguish different 
high-z galaxy populations based on their molecular gas 
content. 
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